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Abstract 


The angles of propagation of the wave fronts 
associated with the duct modes are derived for a 
cylindrical duct with a uniform steady flow. These 
are the angles wliich the normal of the local wave 
front makes with the coordinate axes. The main em- 
phasis is upon the propagation angle with respect to 
the duct axis and its relation to the far* field 
acoustic radiation pattern. When the steady flow 
Hach number is accounted for in the duct, the propa- 
gation angle in the duct is shown to be coincident 
with the angle of the principal lobe of far*field 
radiation obtained using the Wiener-Hopf technique. 
Different Mach numbers are allowed within the duct 
and in the external field. Some interesting results 
of the analysis have implications regarding static 
noise tests and external flow convective effects. 

For static tests with a steady flew in an inlet but 
with no external Mach number the far-field radiation 
pattern is shifted considerably toward the inlet 
axis when compared to zero Mach number radiation 
theory. As the external Mach number is increased 
the noise radiation pattern is shifted away from the 
inlet axis. The theory is developed using approxi- 
mations for sound propagation in circular ducts. An 
exact analysis using Hankel function solutions for 
the zero Mach number case is given to provide a 
check of the simpler approximate theory. 

Introduction 


The angle which the wave front normal vector 
makes with the duct axis is significant since this 
angle governs the location of the principal lobe of 
the far-field radiation. This has been shown in 
semi-infinite rectangular ducts with uniform flow by 
Wright^^' and Cerdclt^) and in circular ducts for 
the no-flow case by Homicz and ^ordi.”' Steady 
flow in the duct will be shown to have an Important 
effect upon the f«r-fielci radiation pattern. 

Also of significance Is the angle which the 
wave front normal makes with the radia'* ^.ordinate 
which is the angle incidence upon tne wall. This 
angle of Incidence Is intimately related to the mode 
cut-off ratio %diich has been shown by to 

correlate the performance of acoustic suppressors. 
The angle of incidence upon the wall is also oecea?. 
•ary for ray tracing techniques as used by Posey 
Ray tracing techniques were also used by Jacques 
to obtain approximate flight effects upon sound ra- 
diating from Jet exLiuat pipes* In refs. 6 and 7 
the connection between the acoustic ray angles and 
the duct modes wer* not considered. It ahould be 
noted when considering trgles of propagation of the 
wave fronts comprising a mocj in s soft-wall duct, 
that these angles will be modified by the acoustic 
liner boundary condition. This effect will be re- 
ported in s separate papei.(^) 

*Head-Acoustlcs Section, Menher AIAA 
**Aerospac^ Engineer, Meiuber AIAA 
tAcoustics Conaultant 


The angles of propagation will be derived from 
the convective wave equation and its solutions. 
Additional convective effects must then be consid- 
ered both in the duct and in the far-field so that 
Inferences can be made regarding the far-field radi- 
ation for both static and external flow conditions. 
These cases are handled by allowing different Mach 
numbers in the duct and in the surrounding medium. 
Approximations will oe used to obtain simplified 
equations for the angles of propagation. Exact 
solutions fer these angles using Hankel function 
soluti >ns are given in appendix A and are used to 
check the approximations made in the main text. 


Symbols 

C modal pressure coefficient, see eqs . (2), 

(A- 2) and (A- 3), N/m^ 

c speed of sound, or the magnitude of the wave 

velocity vector normal to a plane wave 
front, m/sec 

c^ components of vector c, i * x, r, 0, tn/sec 
resultant velocity vector, m/sec 

H^^^ Hankel function of first kind and order ra, 
see eq. (A-4) 

( 2 ) 

Hm Hankel function of second kind and order m, 
see eq. (A-4) 

h magnitude of 

m m 

J Bessel function of the first kind of order m 

- -1 

K magnitude of K, m 

K’ dimensionless local wavenumber, see eq. (A-18) 
K local wave number vector, 

k wave number, CiVc, m"^ 

radial wave number, m“^ 

k combined radial-circumferential wave number 

(a/r^), m-1 

axial wave number, 

k transverse wave number In rectangular duct, 

kg circumferential wave number, m“^ 

M Mach number 

hQ Mach number in duct 

Mach number in surrounding medium 

m spinning mode lobe number 

2 

P scoustic pressure, M/m 

far-field pressure for static tests, N/m^ 
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where s Is the nuisber of lobes for the spinning 
node (circumferential order), \i the radial order 
number of the mode, the axial wave number, 

U) is the circular frequency, Is the outer wall 


radius. 


Is the mode eigenvalue, and is the 


Bessel function of the first kind and order m. 


2 

far-field pressure with external flow, N/m 

outward traveling pressure wave in a cylin- 
drical duct, see eq. (A-2), N/m^ 

inward traveling pressure wave in a cylin- 
drical duct, see eq. (A- 3), N/in^ 

radial coordinate, m 
duct radius, m 

circumferential arc length, m 
time, sec 

axial group velocity, Wsec 

phase velociwy normal to wave front in cylin- 
drical duct, m/aec 

axial coordinate, m 

Bessel function of the second kind and order 
ID 

rectangular coordinate, m 

hardvall duct mode eigenvalue 

circumferential coordinate, radians 

mode cut-off ratio, see eq. (11) 

phase of radians 

m 

wave front propagation angle relative to 
radial coordinate, angle of incidence on 
the wall, deg 

wave front propagation angle with respect to 
circur.fereri ial arc direction, same as 
deg ® 

wave front propagation angle relative to 
axial coordinate, deg 

wave front propagation angle measured from y 
coordinate, deg 

wave front propagation angle relative to 
cir umferential coordinate, deg 


For brevity, the s|i subscripts will be de- 
leted, and it will be understood that a a ingle (but 
quite general) mode is being considered. When 
eq. (2) la inserted into eq. (1) the result Is, 

2 

(^) + “2 - . 2Mj,kk, + I^k2 - (k - 

( 3 ) 

The first term In cq. (3) is a combined radial- 
transverse wave nuinber often denoted by. 


The approximation will now be made that the 
wave fronts behave locally as plane waves propa- 
gating skewed to the coordinate angles of the duct. 
The final approximate equations will be checked by 
an exact but more cumbersome solution using Hankel 
functions in appendix A. The approxiir.ate solutions 
are found to be sufficiently accurate near tiie cyl- 
indrical duct wall where the angle of incidence 
would be used and where moat of the acoustic inten- 
sity usually exists. The i^^ propagation angle can 
be defined by, 
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This combination of wave numbers ia, as will be 
shown later, the cause of the trouble in defining 
some of the propagation angles in cylindrical ducts. 


X constant phase value, see eq. (A- 10), radians 

t resultant axial propagation angle in duct, 

deg 

angle between duct axis and peak of the 
principal lobe of radiation in the fcr- 
field, deg 

CO circular frequency, radians/sec 

Subscripts 

designates quantity for m^ circumferential 
and ^ 1 *-^ radial mode 


i 

cos * - ■ " (5) 

/p! 

which is derived in appendix B, and where i r, 

0, or X and q. Is the angle between the normal 
to the wave front and the i*^^ coordinate axis. 

When there is no flow, eqs. (3) and (4) can be shown 
to yield the usual result. 



k, for • 0 


( 6 ) 


Development of the Propagation Angles 


and then eq. (3) becomec the standard direction 
cosine expression, 


The wave equation in a circular duct with a 
steady flow can be expressed as. 


(1 - 4 ) 
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1 d^p 1 
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where x, r, and 9 are the usual cylindrical 
coordinates, t is ^ Is Mach ntmber, c the 

speed of sound, and V la the acoustic pressure. 
The solution to eq. (..) for the pressure is, 


cos cp^ - for Mjj « 0 (7) 

For finite Mach number however, from eqs. (3) 
and (4). 

( 8 ) 


Axial Propagation Angle 

The angle of propagation with respect to the 
x-axla can be found without further approximation. 
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The axial wave number can be found by rearranging 
eq. (3) to yield 

where the plus sign has been selected for the radi- 
cal to provide the expression for the wave traveling 
in the positive x direction. Eqs. (5), (8) and 
(9) then yield, 

‘f - ^ f) 


cos 


kr 


- 4 


which is consistent with the derivation of Sofrin 
and McCann, but where cut-off occurs when & 1. 

Eq. (10) then becomes, 


cos « 


-Mp + - 1/i^ 


1 - yT - 1/|2 


Note that when « 0 


cos 




1/r 


sin *= 1/1, for Mp « 0 


where 

t V 

* - 

e r 

It is reasonable to assume that 


k|g « k2 + k| 


(16) 


(17) 


which is the result that would be obtained using a 
rectangular approximation to a thin annulus and is 
equivalent to the approach suggested by Cumpsty.'^^' 
Thus using eqs. (4), (16) and (17). 


‘•■iirw 


(18) 


( 10 ) 

where all of the quantities in eq . (10) are known. 

It is convenient here to Introduce the mode cut-off 
ratio, 


and then using eqs. (5) and (8) the other propaga- 
tion angles can be calculated as 


(11) 




and 




kr 


(l - M. "-{ 


(19) 


( 20 ) 


(12) 


Notice that {pj. and 5?g vary with duct radius 
while does not. Introducing the oode cut-off 

ratio (eq. (11)) and letting r « r^, the angle of 
Incidence at the wall can be written as. 


(13) 


(14) 
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( 21 ) 


which is the same relationship as given for the 
angle of the peak of the principal lobe of far field 
radiation as given in refs. 4 and 10. This corre- 
spondence between propagation angle with respect to 
the duct axis and far-field radiation has been noted 
previously for zero Mach number in refs. 1 to 3. 

This correspondence will be shown, in the next sec- 
tion, to hold for cylindrical ducts with the same 
flow in the duct and in the surrounding roediutc. 

Radial and Circumferential Propagation Angles 

As noted earlier, eq. (12) can b: obtained 
without additional approximations. However, the 
angle of incidence on the wall (cp^) and the angle 
to the transverse or circumferential direction (cp^) 
cannot be obtained exactly using the present me^od 
because as shown In eqs. (3) and (4), the radial 
and transverse wave numbers are combined through 
the mode eigen value (a). The transverse wave 
number can be obtained from eq, (2) as follows, 

me - k^re (15) 


It is now easy to see why the cut-off ratio approx- 
imately correlated tlie optimum wall impedance values 
for acoustic liners in ref. 4. For a given Mach 
number and for small m/ci the angle of incidence 
upon the wall is a function only of cut-off ratio. 
The scatter in the above mentioned correlation oc- 
curred mainly for the high circumferential lobe 
number, low radial order modes where m/a was not 
small. Eq. (21) shows this ratio to be an addition- 
al var able when it is not small. Perhaps a better 
correlating parameter for optinum Impedance would 
be the angle of incidence upon the wall rather than 
the cut-off ratio. This has been investigated and 
will be reported in ref, 8. It should be noted that 
the a.-'ft wall boundary condition will alter the 
angle of Incidence from that given in eq. (21). 

Par-Field Radiation Considerations 

As mentioned previously the axial angle of 
propagation has been ahovm to correspond to the 
angle of maximum noise propagation In the far-fleld 
for aero Mach number. However, for the case of 
uniform Mach mimber both in the duct and in the 
surrounding medium, an expression can be written 
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from the results of refs. 2, 


3» end 12 ss. 


co« tp 



( 22 ) 


numerator of cq. (24) 


c„ + cl^ . 


can be shown to be. 



( 25 ) 


where approximate angle for the peak of 

the principal lobe of radiation in the far field. 

The sero Mach number case of eq. (22) has been shown 
to be a close approximation for the location of the 
peak of the principal lobe of the far-iieXd radla* 
tion pattern for unflanged ducts in ref. 3. Saule 

has shcnni that using the duct eigenvalue to 
obtain estimates of the principal lobe peak angle 
can result in a small error for the case of a 
flanged duct radiation pattern and he provides a 
method for applying a correction to this angle if 
it is considered necessary. 

The first obj* ctive here is to show that the 
resultant axisl an Le of propagation in the duct 
with Mach number coincides precisely with the re* 
suit given by eq. (22). This will provide confi* 
dence that the convective effects upon radiation 
are being properly modeled. Then the more inter- 
esting case of different Mach numbers within and 
outside the duct will be derived. 


Resultant Axial Propagation Angle in the Duct 

The geometry of the wave propagation vectors 
is shown in fig. 1, and the derivation which follows 
is considered to be done at the outer wall. Since 
the axial propagation angle does not vary with 

radius in the present approximate analysis it is 
net crucial with regard to the radial position used. 
However in appendix A, the exact solution using 
Hankel functions show to actually be a function 

of radius so the approximate analysis is truly 
valid only near the outer wall. 

In fig. 1, the vector c (speed of sound) is 
normal to the wavefront. The propagation angles 
q) , cp , and Cp^ are the angles between the vector c 
and ti^e coordinate axes and are calculated as in 
eqs. (12), (19), and (20). The components of c 
are calculated from. 


and C|^ is , 



(26) 


Thus eq. (24) becomes 



(27) 


A comparison of eqs. (22) and (27) shows that they 
are identical. The subscript appears on Mp in 
eq. (27) to denote that it is the Mach number in the 
duct. As the wave front leaves the duct, the angle 
will not be changed (refraction effects ne- 
glected) and if the Mach number outside the duct is 
the same as in the duct (as assumed for derivation 
of eq. (22) in ref. 3) then the drift velocity ef- 
fect is also the same as in the ducc. Thus it is 
now evident that for the case which can be checked 
by exact radiation analysis that the resultant axial 
propagation angle in the duct is identical to the 
angle of the peak of the principal lobe in the far 
field. 


Some other Interesting observations can be made 
from the above relationships. The axial component 
of the resultant velocity vector Cj^, as given by 
eq. (25), is the axial group velocity which could be 
derived from eq. (3) by applying 


V • 

gx 



(28) 


• c cos i - X, r, 0 (23) 

The resultant direction of propagation will not be 
normal to the wave front but will be in the direc- 
tion of the resultant vector C|^. The steady flow 
velocity (c Mj)) must be added to Cy^ in a mariner 
similar to that done in refs. 6 and 7, The resul- 
tant vector Cj^ contains the combined effect, due 
to duct Mach number, of the change in wave front 
directions and the drift velocity. 

The angle which is now of interest for compari- 
son with the far-field radiation principal lobe 
peak (eq. (22)) is the angle t between the resul- 
tant vector and the duct axia (x). This angle 
can be defined by, 


cos 



(24) 


The vectors and c^ are considered constant 

In the derivation. After some manipulation the 


Also at mode cut-off (| • 1) the axial group veloc- 
ity is zero and ^ * 90°, In contrast, from 
eq, (12), the axial propagation angle cut- 

off is, 


cos (29) 

As an example for an inlet with Mjj » -0.4, • 

66.4® while for an exhaust duct with Mjj » 40.4, 

■ 113.6. Thus, for an Inlet, the wave front is 
tilted toward the axis %#hlch has Important implica- 
tions upon static test radiation patterns as will 
become more evident in the next section. For the 
exhaust case the wave appears to be going in the 
wrong direction at mode cut-off. It is the angle t 
and not cp^ %ihich governs the sound propagation in 
the duct. This convective drift velocity effect has 
been used in refs. 6 and 7, however, it was not re- 
lated to modal properties. 
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Effect of External Mach Number on 


Far-Field Radiation 


Previous far field radiation theories applica- 
ble to engine inlets have dealt with the 
special cases of zero Mach number or uniform Mach 
number Inside and outside the duct. For an exhaust 
duct, differences in duct Mach number and external 
Mach number across a discontinuous slip layer have 
been handled In an exact analysis by Savkar^^^^ and 
in an anproxlmate ray tracing approach by Jac- 
ques. However, for an engine inlet there are no 
sharp discontinuities along cylindrical surfaces out 
in front of the inlet which would be required for 
these latter two studies to be valid. The following 
analysis is proposed to account for changes in the 
far-field radiation pattern occurring In static 
tests and in tests with external flow. 

In the previous sections the effects of duct 
Mach number upon the angles of propagation in the 
duct have been derived and have been shown to be 
considerable. The propagation angle which has been 
shown to be controlling for the far-field radiation 
(when properly corrected for convective drift veloc- 
ity) is that between the wave front normal and the 
duct axis. As the wave front passes out of the duct 
inlet, only diffraction and refraction cau change 
the direction of this wave front. Diffraction 
scatters the sound and causes the lobed pattern in 
the far-field, but the bulk of the acoustic power 
still radiates In the predesignated direction de- 
termined by the resultant axial propagation angle. 
Refraction through velocity or temperature gradients 
can bend the wave fronts and could possibly be in- 
cluded prior to the application of the following 
analysis. For the present, however, refraction ef- 
fects will be neglected. 

It is thus assumed that the axial propagation 
angle derived in the duct cnviroamenc is also 

valid in the far-field. It Is thus only necessary 
to apply the external drift velocity correction 
which in general will be different from that of the 
duct. 

In th^ far-field, the normal to the wave front 
lies essentially in the radial-axial plane and only 
axial and radial velocity components need to be con- 
sidered in a cylindrical coordinate system. Fig. 2 
shows a sketch of the velocity vectors. All that 
needs to be done is to apply the external field 
(flight or static) drift velocity correction to the 
vector c. Recall that ' 4 :^ is unchanged from its 
value in the duct and is given by, 



The principal lobe of radiation can be expected at 
an angle where. 


cos tp 


«iR 


( 31 ) 


After some manipulations this angle can be expreased 

by, 



( 32 ) 


If the external velocity is the same as the duct 
steady flow velocity, then eq. (32) reduces to the 
•pedal case of eq. (27). Eq. (32) can be consid- 
ered as an equation describing the effect of a wind 
tunnel velocity since In general Mpo Mp. 


Far-Field Angles for Static Tests 

An important special case of eq. (32) Is that 
of Moo * 0 which is applicable fox static engine 
terts. Eq. (32) then reduces to 



which is of course the same as the relationship 

given by eq. (12). Sample calculations using 
eq. (33) are shown in fig. 3. Of particular inter- 
est are the modes near cut-off. If the duct Mach 
number is zero, the principal lobe of far-field ra- 
diation will occur at 90*^. Ibis of course is also 
predicted by previous radiation theories when Mj) * 

0 or when Mp * However, as the inlet duct 

Mach number increases, the near cut-off modes propa- 
gate more toward the axis. For example at a typical 
inlfjt Mach number of -0.4, * 66.4°. At Mp * 

-0.8, % “ 36.7° which must account lor at least 
soDie of the apparent sideline attenuation of a near 
soJiic inlet static test. For modes above cut-off 
the same trends are observed with a shift of the 
sound radiation toward the axis. 

The results shown in fig. 3 are qualitatively 
corroborated by acoustic suppression results for 
the blade passage frequency as reported in refs. 15 
and 16, ^The maximum suppression occurred at 70 
(only 10 intervals measured) from the inlet at full 
engine speed with a Mach number Mp * -0.375. Fox 
these static tests the blade passage frequency 
should be rich in modal content near cut-off. Pre- 
vious radiation theories would thus predict a peak 
attenuation near 90 for these most easily atten- 
uated modes near cut-eff. The present radiation 
theory, eq. (33) or fig. 3, would predict that peak 
attenuation should occur near 68^ which is in good 
agreement with expetlmental results. Also apparent 
in the data of refs. 15 and 16, is the fact that as 
engine operating speed and the Mach number is re- 
duced the peak attenuation moves more toward the 
sideline which la also in agreement with the pres- 
ent theory. 


Far-Field Radiation with External Flow 

Fig. 4 ahows the far-field radiation peak for 
near cut-off modes when the aurrounding medium is 
also moving which would simulate a wind tunnel test. 
(3’’' as used ';lth • 1 for various duct Mach 
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niubbers (Hp) and external flow Mach number (%>)- 
When Moo 0 ^he results are the same as in fig. 3. 
However, as the tunnel speed increases (1^ pre- 
sumed to remain nearly constant) the sound radia- 
tion moves back toward the sideline. For example, 
sample external velocities irinxilating approach and 
takeoff conditions are shown on fig, 4. For the 
approach condition • *0.4) a atatic test would 
show the near cutoff modes propagating at about 
66° while in a wind tunnel this radiation would oc- 
cur near 78° for hi* - -0.2. The difference be- 
tween a take-off static and tunnel test would be 
15° (from 53° to 68°;. 

Since there is a difference between static and 
wind tunnel far-field radiation patterns even for a 
specific mode, a correction should be applied to 
wind tunnel data before it Is used in a fly-over 
calculation. Assuming that an inflow control de- 
vice has been used in the static tests so that 
flight modal structure has been simulated) the dif- 
ference between the static test angles (cp ) and the 
tunnel Lest angles (>'p) must be reconciled. An ex- 
pression can be derived relating these two angles 
by eliminating the cut-off ratio (|) between eqs . 
eqs . (12) and (32) or more simply by considering 
the geometric relations shown in fig. 2, The final 
expression is, 


(M^* + cos CJ^) 
y 1 + m 2 + 2M„ cos ct. 


(3A) 


This is equivalent to the expression of ref. 17 
which deals with external convective effects upon 
noise radiation. This expression shows that those 
modes which propagate at far-field angle in a 

static test, will propagate at angle in a wind 

tunnel test. This correction can be considerable. 
For example, with a tunnel speed of only Ftp “ -0.2, 
the static data measured at 50°, 60°, and 70° should 
be shifted to 60^^, 71°, and 81.5® for a wind tunnel 
test. This will have the effect of shifting the 
usually higher sound pressure level data of these 
foi-ward angles back more toward the sideline. It 
was assumed in the above discussion that the princi- 
pal lobes dominate the radiation pattern and that 
even for a muitimcdal pattern the shift will thus 
be properly described by eq. (34). 

Technically a correction should also be made 
for the sound pressure level when the angle shift of 
eq. (34) is made. By maintaining acoustic power 
and accounting for the principal lobe width change 
(solid angle change) which occurs with the angle 
shift the following expression can be derived, 


(1 + + 2M„ cos tp 

00 QO 00 * 

— es - " ■ ■ 

(1 + fL cos tP,) 

O "X 


(35) 


For tunnel Mach numbers between -0.2 and -0.4 the 
sound pressure level shift for the static test 
angles between 50° and 90° is less than 1.5 dB. 

Thu/i for the sideline angles of most interest the 
sound pressure level correction Is probably not 
si^^nif leant, although the angle ahift of eq. (34) 
should definitely be considered. It should again 
be noted that refraction effects within the velocity 
grsdients near the Inlet for static vests have not 
been considered here. These effects should not be 
great, however, for the near sideline angles of 
most Interest. 


While the radiation theory presented In this 
paper appears to receive preliminary confirmation 
from static engine tests, there is not yet any con- 
firmation of the wind tunnel effects since there 
were no dsts svsilable. Careful testing would be 
required to insure that the results were not masked 
by other factors. For example, the modal structure 
of a fan source would be expected to be different 
for static testing than fer wind tunnel testing , ' 
Inflow control devices (screens) would be required 
to provide a static test noise source which would 
have a chance to simulate the modal content produced 
by s turbofan In a wind tunnel. 

It should be noted that the corrections dis- 
cussed above must be made for data obtained in a 
wind tunnel test since it is the directivity pattern 
described as a function of and not which 

is projected to the stationary observer in a simu- 
lated fly-over calculation. 

Concluding Remarks 

The angles of propagation for the wave front 
making up a duct mode have been presented here with 
Mach number in the duct. Approximate equations have 
been derived to provide simple utilitarian expres- 
sions. These expressions are valid only near the 
outer wall which is the most important region since 
the bulk of the acoustic intensity is located there 
and this is also where incidence angles would be of 
interest. Exact solutions using Hankel functions 
are given in appendix A and these corroborate the 
approximate solution accuracy near the outer wall. 

The main emphasis of this paper was to use the axial 
propagation angle to infer information about the 
far-field radiation pattern. The resultant axial 
angle of propagation in the duct was shown to agree 
exactly with the peak of the principal lobe of far- 
field radiation obtained from formal radiation 
calculations when the Mach number is uniform everv- 
where. The present solution was then extended to 
cover the case of different Mach numbers inside and 
outside the duct for which exact calculations have 
not been available for engine inlet configurations. 
The new analysis shows that for static engine tests 
the inlet radiation can be expected to be shifted 
considerably toward the inlet axis over that obtain- 
ed from previous analysis. This static test radia- 
tion shift can be shown to have preliminary verifi- 
cation. An external Mach number convective effect 
is also predicted which would shift the sound radia- 
tlcxi back toward the sideline (compared to a static 
test) as the external flow velocity increases. The 
refraction of the sound by the flow gradients near 
the inlet could possibly be Included in the analy- 
sis, but this effect has not been included here. 

Appendix A 

Cylindrical Wave Synthesis of a 
Ducted Spinning Mode 

The exterior radiation field of an Infinite 
cylinder sustaining standing circumferential surface 
oscillations, synchronous in the axial direction, 
is shown by Morse and Ingard^^^' to be described by 
the Hankel function combination, cos(m 6) Hj^^^(kr) 
exp (-itDt). As a separate matter in the same ref- 
ference It la demonstrated that the modal compo- 
nents of the pressure field inside a hardwall rect- 
angular duct may be eynthesized by complementary 
pairs of inward and outward plane wave trains. 
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These concepts may be combined and extended to 
describe the field of a spinning mode in a cylin* 
drical duct. 

Let the (complex) pressure - Pn,(x,©,r,t) 
in the duct be expressed as the sum of a pair of 
cylindrical vaves: 


P - P^ + P^ 
tn m m 


The radially outward wave is given by 


pi . 1 c 
m 2 aw m 


l(k^x+m^-:Dt) 


and the matching Inward wave is 

•> 1 i(k *+m9-ajt) 


2 au m \tJ 


The eigenvalue a as in the main text should con- 
tain the subscripts m,u hut these have been 
dropped for brevity. and H^2) are Hankel 

functions of the first and second kinds: 


+ iY 

m m m 


Inward and outward wave designations are apolied as 
a consequence of the far- field nature of and 

which is disclosed by the asymptotic behavior 
of the Bessel functions; 


n — ifk x+in0:t~-OD 

|r mr 


With these substitutions, the general outward wave 
expression, eq. (A-2) becomes 


ifk x4me-Hi Va.q 

r\ L * “VV J 


i C h (^)e 
m 2 inA ®Up/ 


The local wave direction and phase velocity are 
obtained by examining the exponent 


X - k x + m6 + 0 (“ 
X a\r 


A plane that Is locally tangent to the wave 
front will be governed by the requirement that the 
phase, X is conserved over a small displacement 
occurring during a corresponding small time inter- 
val. For this phase conservation dX * 0 or 


dx * k dx + mde + t- dr - cijdt * 0 

X dr m\r i 

\ 0 / 


Writing 


(0 ‘ 


and introducing the circumferential arc length coor- 
dinate 

S - re (A- 13) 

the equation of the local phase plane becomes 

k^dx * 2 ds + ^ ‘ 


(A-lAa) 


Two comments may be helpful. First, it is ob- 
vious from eqs. (A-2), (A-3), that the sum of ^s^^ ^r^^ (A- 14a) 

and gives the duct mode field of eq . (2). 

Second, the opposite signs applicable to re- Xhe coefficients of the differential coordinates 

quire the wave to move radially inward in order to in (A-14) are the axial, tangential, and radial 


conserve phase. 

The near- field behavior of the outward wave 
may be found by representing the Hankel function 


may be i< 
Hiih in 


polar form: 




where the amplitude 




and ^ , the phase of H , is 

ID ID 


t.n - 


coinponents_^of the local wave number vector, de- 
noted by K. On dividing (A-14) by the magnitude 
of K, denoted by K there follows: 


cos dx + cos ds + cos dr 


cos 

Cp B 

k /K 


x 

x 

cos 

m 

■ 

cos t 

cos 


— 9> 


r 

r t 
o 
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£q. (A<-I5) Is the normal form of a plane with 
direction angles :pg, traveling at a local 
phase velocity Vpj, * a>/K. If K is expressed as: 



the magnitude of the phase velocity becomes 


Vph 


K 


CJD 

kK‘* 


c 

K' 


(A- 19) 


Thus the phase velocity differs from the speed 
of sound by the factor I/K' . It can be shown that 
K’ > 1 and approarnes 1 as (ar/r^) », so that the 

near*fleld phase velocity is always subsonic. This 
means that the sum of the squares of the wave num- 
bers kjj, kg, kr, add up to , which is larger 
than the square of the ordinary wave number, k ■ 
Cd/c. 

Further work requires evaluation of 
the derivative of the phase of the Henkel function 
H^^)(ar/ro). It is found that 



This expression allows compulation of the l*^cal 
wave speed and direction angles for a wave at any 
radius usinc eqs. (A-16) to (A- 19), 

some ralculstlons were made for the wal' radius 
location r » Tq, in the Mj) * 0 case, for cutoff 
ratios of 1 and 2. The tabulation in table Al pre- 
sents values of and the wall incidence angle 

cpj. (eq. (A-16)), For comparison, values of 


computed from the main text eq. (21), which is 
based on Che approximation K • k or K* » 1 i are 
Included. These values show that the greatest de- 
partures from ’’conventional" or far- field behavior 
occur near cutoff for the lowest order m and p 
modes. The departures, however, are quite minor: 

For the "worst" case, the (1,0) mode st cutoff, the 
phase velocity departs from sonic by only 6 percent 
and the exact anU approximate wall incidence angles 
differ by about 2^. Thus the essential assumption 
used in deriving the results of the main text, 
namely that for * 0 the squares of the axial, 
tangential, and radial wave-number components sum 
to the square of the ordinary wave number, k « U)/c, 
at the wall radius, is in excellent practical agree- 
ment with results of the exact analysis. This as- 
sumption implies that the eigenvalues, which are 
independent of are related by kg + k|y • 

k|e - (a/ro)2. 

Finally, it may be of some interest to note 
that the situation changes radically as the duct 
axis is approached more closely. Calculations show 
the following (m ^ 0): 


K' - « 


cp - 90° 

X 

'tj. - 90° 


OF 










v.v. 

cio: 


Thus iiear the axis the local wave velocity is very 
small and is essentially circumferential in direc- 
tion, This contrasts with the wall behavior where 
the phase velocity is almost sonic and the wave di- 
rection is predominantly axial-radial. 


TABLE Al, - VAUIES OF K’ AND FOR « 0, r » r^ 


Mode 

(niiP) 

C 

utof f ratio | 

« 1 

Cutoff ratio i 

* 2 








K- 

::;^(exact) 

'4^^(approx) 

K’ 

:;^(exact) 

q.^(approx) 



deg j 


deg 

1.0 

1.0633 

30.7 

32.9 

1.0162 

63.3 

65.2 

1,1 

1.0050 

10.8 

10.8 

1.0013 

60.5 

60.6 

1,2 

1.0018 

6,7 

6.7 

1.0005 

60.2 

60.2 

2,0 

1.0462 

38.8 

40.9 

1.0117 

66.2 

67.8 

2,1 

1.0043 

17 .3 

17.4 

1.0011 

61.4 

61.5 

2,2 

1.0019 

11.6 

11.6 

1.0005 

60.6 

60.7 

4,0 

1.0334 

46.7 

48.8 

1.0084 

69.4 

70.8 

4,1 

1.0030 

25.4 

25.8 

1.0008 

63.1 

63.2 

4.2 

1.0012 

18.4 

18.4 

1.0003 

61.6 

61.7 

8,0 

1.0232 ' 

54.1 

56.0 

1.0059 

72.7 

73.8 

8,1 

1.0030 

34.4 

34.5 

1.0002 

65.6 

65.7 

8,2 

1.0010 

26.7 

26.7 

1.0002 

63.5 

63.5 

16,0 

1.0158 

60.7 

62.4 1 

1.0040 

75.7 

76.6 

16,0 

1.0018 

43,4 

43.5 

1.0004 

68.7 

68.7 

16,2 

1.0010 

35.8 

1 

35.8 

1.0002 

66.1 

66.1 
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Appendix E 

Wave* front Propagation Angles vlth Flow 

Itiis appendix is Intended to provide a siiBpIi- 
fied example of obtaining the angles of propagation 
of the wave fronts in a duct with flow. To serve 
this purpose a simple two dimensional rectangular 
geometry will be used. The wave equation solution 
can be given as, 

l( 0 )t-k X) 

p • e cos kyY 


• >PJGINAL PAGE IS 
I' POOR OIJAEITY 

and some rearrangement of eq. (B* 6 )' yields, 

+ ^1 - (1 - mJ ) 


-Md 


(B-8) 


The wave number ky is not a function of lip, but 
alnce the denominators of eqs. (fi-4) and (B*5) are 
functions of h%), then both (p^ and depend 

upon duct Mach number. 
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Figure 4. • Effect of external and duct velocities on far- 
field propagation angle for near cutoff modes, 


-N0Rr/,AL TO WAVE 


WAVE FRONT- 


u/ 

■ , . 






Figure 2 , - Fcr-field propagation angles with flow. 


Figure Bl. - Sketch cf propagation angles and wavefront slope m a two 
dimension rectangular duct. 
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Figure 3. - Convective effect upon tar-field radiition 
principle lobe angle for inlets under static test 
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